ONR-URI  Composites  Program 
Technical  Report  No.  91-02 


UIUC-NCCMR-91-02 


AD-A234  102 


A  MODEL  FOR  WET  FILAMENT  WINDING  OF  COMPOSITE 
CYLINDERS  WITH  ARBITRARY  THICKNESS 

A.  Agah-Tehrani*  and  Hong  Teng** 


January,  1991 


National  Center  for  Composite  Materials  Research 
at  University  of  Illinois,  Urbana  -  Champaign 
A  DoD  University  Research  Initiatives  Center  funded  by  the 
Office  of  Naval  Research,  Arlington,  VA 


Assoc.  Mem.  ASME 
Graduate  Assistant 


0\\  ^  \bU> 


A  Model  for  Wet  Filament  Winding  of  Composite  Cylinders  with 

Arbitrary  Thickness 

A.  Agah-Tehrani1  and  H.  Teng2 

Department  of  Theoretical  and  Applied  Mechanics 
University  of  Illinois  at  Urbana-Champaign 
Urbana,  IL  61801 

ABSTRACT 

A  continuum  consolidation  model  is  purposed  for  analyzing  the  filament  winding  of  thermoset 
composite  cylinders  with  arbitrary  thickness.  The  model  takes  account  of  the  variation  of  both  the 
instantaneous  stiffness  and  the  permeability  of  the  mixture  with  fiber  compaction.  Due  to  resin 
filtration,  the  structure  of  the  resulting  initial-boundary  value  problem  is  similar  to  that  of  a  moving 
boundary  problem.  Based  on  this  analogy,  a  finite  difference  scheme  is  devised  for  the  solution  of 
the  problem.  For  the  case  of  winding  onto  a  rigid  mandrel,  the  results  point  to  the  existence  of  an 
active  and  a  passive  zone  of  consolidation.  The  results  further  indicate  the  possibility  that  during 
hoop  winding  of  relatively  thick  cylinders,  the  tension  can  be  completely  lost  in  the  portion  of  the 
passive  zone  away  from  the  mandrel. 
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INTRODUCTION 

The  filament  winding  is  a  manufacturing  process  which  can  be  used  for  forming  composite 
cylinders  with  high  strength.  The  essence  of  this  processing  method  is  to  successively  add  layers 
with  different  winding  angles  onto  a  moderately  stiff  cylinder  which  is  invariably  circular  (mandrel). 
The  layers  may  be  wrapped  in  adjacent  bands  or  in  repeating  patterns  which  eventually  cover  the 
surface  of  the  mandrel.  The  winding  angle  may  vary  from  longitudinal  (parallel  to  the  axis  of  the 
mandrel)  to  circumferential  (  perpendicular  to  the  axis  of  the  mandrel).  The  binder  is  a  thermoset 
resin  which  can  be  epoxy,  polyester,  or  vinyl  ester.  If  the  thermoset  resin  is  mixed  with  the  fiber 
bundle  during  the  winding  stage,  one  has  wet  winding.  Alternatively,  if  pre-impregnated  B-staged 
form  of  the  mixture  of  resin  and  fiber  bundle  is  used,  one  has  dry  winding  (  for  an  exhaustive 
discussion  of  filament  winding  see  Shilbey,  1982).  For  a  given  composite  system,  the  primary 
control  variables  during  winding  can  be  identified  as:  the  tension  applied  to  each  layer,  the  velocity 
of  the  machine  which  traverses  along  the  axis  of  the  mandrel  (cross-head  velocity),  and  the  angular 
velocity  of  the  mandrel  (Calius,  Springer,  1985). 

Following  the  winding  stage,  the  wound  cylinder  is  initially  cured  at  elevated  temperatures 
so  as  to  cause  cross-linking  of  the  resin,  and  then  it  is  cooled  down  to  ambient  temperature. 
Alternatively  one  can  envisage  the  case  where  curing  is  accomplished  during  winding.  In  such 
case,  the  rate  of  heating  the  mandrel  should  also  be  considered  as  one  of  the  control  variables.  In 
order  to  predict  the  occurrence  of  damage  in  the  final  product  (which  is  often  in  the  form  of  the 
fiber  waviness  or  the  transverse  cracks  in  between  the  layers),  there  is  a  need  for  a  i enable  method 
for  predicting  the  state  of  the  composite  cylinder  at  the  end  of  each  of  the  three  stages  of  the 
manufacturing  process:  1  )winding,  2)curing,  and  3)final  cool-down.  For  both  dry  and  wet  winding, 
this  requires  prediction  of  the  stresses  which  are  generated  during  the  winding  process.  During  wet 
winding,  one  must  also  determine  the  radial  distribution  of  the  fiber  volume  fraction  which  arises 
as  the  result  of  the  settling  of  fiber  bundles  in  the  liquid  resin. 

By  assuming  that  successive  layers  are  wrapped  simultaneously  on  the  mandrel,  one  can 
utilize  the  existing  models  for  stress  analysis  of  tape  winding  ( Willett  and  Poesch,  1988;  Lin  and 
Westmann,  1989)  in  order  ro  simulate  the  dry  winding.  In  so  doing,  one  is  neglecting  the  fact  that 
in  reality  the  layers  are  wrapped  sequentially  next  to  each  other  resulting  in  the  axial  variation  of 
the  outer  radius.  The  analogy  can  be  further  completed  by  observing  that  dry  winding  is  a  generalized 
plane  strain  problem  while  tape  winding  is  a  plane  stress  problem. 

Irrespective  of  the  thickness  of  the  wound  cylinder,  theoretical  modelling  of  wet  winding 
must  take  account  of  the  resin  filtration  through  tne  fiber  bundles  which  is  induced  by  the  existing 
pressure  gradient.  Recently  Calius  and  Springer  (1990)  have  outlined  a  model  for  simulating  the 
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entire  stages  of  wet  winding.  Their  submodel  related  to  the  winding  stage,  which  neglects  radial 
resistance  of  the  fiber  network,  is  based  on  considering  the  force  equilibrium  and  the  conservation 
of  mass  within  each  layer.  As  it  was  duly  noted,  this  model  is  only  appropriate  for  thin  cylinders 
where  the  increase  in  the  radial  stiffness  due  to  the  consolidation  of  the  fiber  bundles  can  be  neglected. 
Within  the  framework  of  a  discretized  consolidation  model,  Bolotin  et  al.  (1980)  have  included  the 
elasticity  of  the  fiber  network  by  introducing  the  term  associated  with  the  radial  pliability  in  the 
equilibrium  of  the  layers  whose  thickness  is  below  the  critical  value  at  which  the  radial  resistance 
arises. 

Any  comprehensive  model  of  wet  winding  must  take  account  of  two  important  aspects  which 
are  associated  with  fiber  compaction:  i)  increase  in  the  radial  stiffness  of  the  composite  (Gutowski, 
et.  al;  1987a, b),  and  ii)  reduction  in  the  permeability  (Gutowski,  et.  al;1987a,b,  Gebart;  1990).  The 
significance  of  these  attributes  become  particularly  apparent  when  the  thickness  of  the  cylinder 
becomes  large.  The  primary  aim  of  this  paper  is  to  provide  a  continuum  consolidation  model  for 
analyzing  the  wet  winding  of  thermoset  composite  cylinders  with  arbitrary  thickness  which 
incorporates  the  above  characteristics.  Such  a  model  is  developed  by  utilizing  the  rate  form  of  the 
Biot’s  consolidation  mixture  theory  (  Biot;  1941).  In  comparison  with  the  existing  discretized 
consolidation  model  purposed  by  Bolotin  et.  al.  (1980),  the  present  formulation  has  two  distinct 
advantages:  One  is  that  it  is  based  on  a  continuous  model,  and  the  other  is  that  it  takes  account  of 
the  reduction  of  permeability  with  fiber  compaction.  Following  the  description  of  the  key  features 
of  the  model,  the  details  of  the  numerical  calculation  and  interpretation  of  the  results  will  be 
presented.  The  article  will  be  concluded  by  a  discussion  on  the  importance  of  including  the 
dependence  of  the  radial  resistance  and  permeability  on  fiber  volume  fraction  in  the  analysis  of  wet 
winding. 

I.  Development  of  the  Model 

The  composite  cylinder  is  assumed  to  be  a  mixture  of  incompressible  fluid  (resin)  and  solid 
(fiber  bundle).  The  fiber  bundles  are  assumed  to  be  in  the  hoop  direction.  The  axial  flow  of  the 
resin  due  to  the  time,  delay  in  winding  along  the  axis  of  the  cylinder  will  be  neglected;  i.e  it  is 
assumed  that  all  layers  along  the  axis  of  cylinder  are  laid  simultaneously.  Moreover,  the  angular 
variation  of  the  outer  radius  due  to  the  helical  nature  of  the  winding  will  also  be  neglected,  see 
Fig.l.  Because  of  these  two  assumptions,  the  only  non-vanishing  component  of  the  velocity  of 
each  phase  will  be  the  radial  velocity,  and  the  problem  becomes  one  dimensional  in  space.  In  order 
to  facilitate  the  analysis,  the  following  non-aimensionai  group  of  parameters  are  introduced 
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r=R,x  ,  us=Rtus 


,  Hj~Rl  uf 


y  = 


(Orr>G*»P)  -CJ&rr,GmP)  » 


(1) 


where  R,  is  the  outer  radius  of  the  mandrel,  us  and  uf  are  radial  displacements  of  the  solid  and  fluid 
phase,  C°r  and  C$r  are  the  moduli  of  the  fiber  network  which  has  reached  the  maximum  achievable 
fiber  volume  fraction,  and  C<^  is  the  hoop  stiffness  of  an  individual  fiber.  arr  and  6^  are,  respectively, 
the  radial  and  the  hoop  stress  components  of  the  mixture,  p  is  the  partial  pressure  of  the  resin,  t 
denotes  time,  and  k°  is  the  permeability  as  the  fiber  bundle  passes  through  the  resin  bath. 


By  utilizing  the  above  relations,  one  can  determine  the  following  expressions  for  the  non- 
dimensional  velocities  of  the  fiber  network  and  the  resin  denoted  respectively  by  ws  and  w/5 


w. 


dus  tf  diis 
9y  /?,  dt 


duf  tfddf 

w^=RtlF 

with  tf  =  RfhrfCm  as  the  characteristic  time  for  filtration. 


(2) 

(3) 


The  dependence  of  the  instantaneous  moduli  of  the  fiber  bundle  C4n  and  the  instantaneous 

permeability  of  the  fiber  network  on  the  fiber  volume  fraction  0,  are  assumed  to  be  of  the  form 
given  below 


(crr,cj  =  F(<y(C,o  ,  Cee^c, 

(4) 

K=K°C  (<!>,)  • 

(5) 

The  values  of  the  moduli  C°  and  Cr°6,  and  the  functions  F(tys)  and  G(< J>,)  must  be  determined  from 
experiments  on  fiber  compaction  (Gutowski,  et.  al;  1987b,  Gebart;  1990).  It  should  be  noted  that 
for  simplicity  we  have  assumed  identical  dependence  of  Crr  and  on  the  fiber  volume  fraction 


In  terms  of  the  apparent  density  of  the  phases  (p,  =  m^/V^  ,  i  =  s,f),  the  conservation  of  mass 
tor  the  fluid  and  the  solid  phase  can  be  written  as  (Bowen;  1976,  Kenyon;  1976) 


3  It  is  important  to  note  that  these  are  the  velocities  of  a  continuum  point  which  contains  the  two 
phases. 
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with  V  denoting  the  gradient  operator. 


By  assuming  that  the  densities  of  both  solid  and  fluid  phases  remain  constant  during  defor¬ 
mation,  (6)  and  (7)  can  be  manipulated  to  obtain  the  following  expressions  for  the  conservation  of 
mass  in  terms  of  the  volume  fractions  of  individual  phases  (<{>,  =  p.V./m, ,  i  =f,s ) 


dy 


+  V 


*  [fjj/Wy]  =  0 


dy 


+  V 


•  M>,wJ  =  o 


(8) 

(9) 


Addition  of  the  above  two  relations,  and  usage  of  the  identity  $/  +  <}),  =  1  leads  to  the  final 

desired  form  of  the  conservation  of  mass  which  involves  only  the  velocity  of  the  solid  and  the 
filtration  flux  q  =  tyf(wf  -  wj  (the  surface  integral  of  q  over  the  outer  radius  of  the  cylinder  is  the 
amount  of  the  fluid  leaving  the  mixture).  By  enforcing  the  condition  that  no  resin  passes  through 
the  mandrel  (<?(l,y)  =  0),  The  resultant  expression  can  be  integrated  to  obtain 

q=-ws+-8(y)  .  dO) 

where  g  (y )  is  the  velocity  of  the  interface  between  the  mandrel  and  the  composite.  We  assume  that 
the  mandrel  undergoes  small  deformation  so  that  changes  in  the  outer  radius  of  the  mandrel  can  be 
neglected. 

An  important  result  of  the  conservation  of  mass  analysis  is  the  ability  to  determine  the  evo¬ 
lution  equation  which  governs  the  change  in  the  fiber  volume  fraction.  To  obtain  such  relation, 
one  can  simply  rearrange  the  terms  in  (8)  to  get 

d§} 

^  =  -<UV.wJ  ,  (11) 

where  d/dy  represents  the  total  derivative.  The  above  equation  by  itself  does  not  contain  the  property 
tiiat  the  fiber  volume  fraction  should  not  exceed  the  achievable  volume  fraction.  In  other  words, 
if  v,'t  <  0,  the  above  expression  does  not  provide  any  upper  bound  on  the  value  of  <J>,.  However,  it 
is  understood  that  although  the  achievable  fiber  volume  fraction  <j>0  varies  with  the  packing  geometry, 


it  will  never  exceed  one.  This  constraint  implies  that  as  ^  — >  <j>a,  ws  — » 0-  As  it  will  be  indicated 
later,  this  constraint  will  be  satisfied  if  the  variation  of  permeability  with  fiber  volume  fraction  is 
taken  into  account. 


Since  during  wet  winding  the  radial  resistance  of  the  layer  just  being  wound  on  the  mandrel 
is  very  small,  the  resin  material  can  flow  through  the  fiber  network  leading  to  increases  in  both  the 
fiber  volume  fraction  and  the  radial  stiffness  of  the  composite.  Because  of  the  finite  changes  in 
the  fiber  volume  fraction,  one  needs  to  distinguish  between  the  initial  and  the  instantaneous  con¬ 
figuration  of  the  mixture.  In  so  doing,  we  can  either  follow  the  motion  of  the  fluid  panicle  (resin) 
or  the  solid  panicle  (fiber  network).  Since  the  final  radius  of  the  composite  cylinder  depends  on 
the  position  of  the  last  fiber  layer  which  is  being  wound  on  the  mandrel,  the  instantaneous  con¬ 
figurations  of  the  mixture  and  the  solid  phase  are  considered  to  be  identical.  Thus,  or,(x,y)  and 
Gm{x ,  >' )  the  stress  components  of  a  mixture  point  at  "time"  y  which  is  occupying  the  position 
x  of  the  solid  phase.  In  conjunction  with  this,  we  adopt  an  "Updated-Lagrangian”  formulation  to 
describe  the  motion  of  the  mixture,  that  is  at  each  instant  of  time  the  current  configuration  of  the 
solid  phase  is  considered  to  be  the  initial  configuration  of  the  mixture. 


Due  to  the  variation  of  the  outer  radius  with  time,  the  constitutive  behavior  of  the  solid  phase 
and  Darcy’s  Law  are,  respectively,  expressed  in  the  following  rate  form 


Grr  +  p  =/r«>,) 


L  dx  x  J 


dw,  w, 

0„+p=F®s)^  +  hj  , 


(12) 

(13) 


where  (  )  =  d/dy,  and  p  is  the  resin  pressure  in  the  current  configuration. 


(14) 


The  presence  of  the  term  G  (<(>,)  in  ( 14)  implies  that  as  G  (<j>f)  — »  0,  the  filtration  flux  approaches 

zero  (  fiber  compaction  comes  to  a  halt).  For  the  case  of  the  rigid  mandrel  (g(y)  =  0),  it  immediately 
follows  that  the  velocity  of  the  solid  will  also  approach  zero.  Thus,  we  have  shown  that  at  least 
when  the  mandrel  is  rigid,  considering  the  variation  of  the  permeability  with  0,  leads  to  a  bound  on 
the  fiber  volume  fraction.  Further  limitation  on  the  increase  in  compaction  can  be  achieved  by 
incorporating  the  variation  of  radial  stiffness  with  fiber  volume  fraction.  This  can  be  easily  explained 
by  observing  that  the  compaction  will  be  suppressed  once  it  has  reached  a  level  where  the  radial 
stiffness  of  the  fiber  network  is  no  longer  negligible. 

Substitution  of  (10,12,13,14)  into  the  rate  equilibrium  equation  in  the  radial  direction 
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05) 


Orr  ^86 
x 


=  0  , 


leads  to  the  following  parabolic  partial  differential  equation  (PDE)  for  wJ 


d_ 

dy 


w. 


G®,)] 


=  F($S)}\ 


cfw,  1  dw, 

dx2  +  x  dx 


ws 

X 


+ 


dF®t)db 

d§,  dx 


f.  dws  tO 

i  a 

g(y) 

h*+p7j 

+  x9y 

Lo(<t>,)J 

(16) 


The  value  of  the  function  F  (<{),)  will  be  very  small  (close  to  zero)  near  the  outer  layer  (where 

fiber  compaction  has  not  taken  place),  while  in  the  inner  layers  where  enough  compaction  has 
occurred  it  will  have  a  finite  value.  Because  of  this,  the  above  equation  is  an  example  of  a  singularly 
perturbed  PDE. 


2.  Boundary  Conditions 

The  formulation  of  the  problem  will  be  completed  by  providing  the  necessary  and  sufficient 
boundary  conditions  for  solving  (16).  The  first  boundary  condition  derived  will  be  that  associated 
with  the  rate  of  change  of  the  fluid  pressure  at  the  instantaneous  outer  radius  xa.  This  can  be 
accomplished  by  first  substituting  (10)  into  (14),  integrating  the  resultant  equation,  and  enforcing 
the  boundary  condition  that  at  p  (jc  =  xOJ  y )  =  0  for  all  the  time  y .  The  end  result  will  be 


1 

G®,) 


w,(4,y) 


(17) 


Based  on  the  above  expression  and  (16),  the  material  derivative  of  resin  pressure  will  be 
w,fo,y)  1  g(y) 


dx0 

Mx,y)=-^ 


L  G(C)  J»C(«)J 


•JC 


r;i  A-pSiv,  <k-F(<j>,)p 
F  «>,)  — - - - w. 


%  a* 


d$  + 


F(k) 


dw.  w, 

x-^+p- 

dx  x  . 


(18) 


where  4>”  is  the  fiber  volume  fraction  in  the  layer  just  being  wound. 


Following  Lin  and  Westmann  (1989),  we  let  the  amount  the  material  added  to  the  cylinder 
during  the  winding  of  each  layer  be 
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d(Vol/unit  length)  =  2nxadx  =  w0{y)dyH  ,  (19) 

with  H  as  the  thickness  of  each  fiber  network  just  being  wound,  and  wc  as  the  speed  of  the  last 
layer.  In  considering  different  possibilities  for  the  variation  of  w0  with  "time"  y,  one  can  recognize 
two  special  cases:  one  in  which  w„{y)  =  const,  and  the  other  in  which  the  angular  velocity  is  held 
constant. 


If  the  amount  of  resin  filtration  is  small,  then  changes  in  the  geometry  can  be  neglected  and 
(19)  can  be  integrated  to  determine  the  outer  radius  as  a  function  of  time.  Indeed,  this  is  the  approach 
adopted  during  the  elastic  or  viscoelastic  analysis  of  winding,  (Willett  and  Poesch,  1988;  Lin  and 
Westmann,  1989).  In  the  current  problem  where  resin  filtration  leads  to  changes  in  the  position  of 
the  fiber  network,  equation  (19)  will  be  interpreted  with  respect  to  instantaneous  current  configu¬ 
ration  of  the  fiber  network.  In  this  case  (19)  leads  to 


dx0  w0H 
dy  2nx0 


(20) 


The  fact  that  ihe  evolution  of  the  outer  radius  of  the  composite  ( the  position  of  outermost 
fiber  network)  is  not  known  a  priori  at  each  moment  in  time  indicates  that  the  wet  winding 
problem  is  in  some  sense  an  example  of  a  moving  boundary  problem,  and  one  needs  to 
sequentially  integrate  the  velocity  equilibrium  equation  (16)  for  each  instant  of  time  in  order 
to  obtain  the  relevant  information  at  the  end  of  the  winding  process.  To  accomplish  this  one 
needs  two  boundary  conditions  on  either  the  velocity  or  the  strain  rate  of  the  solid,  or  a  mixed 
combination  of  both  (recall  that  (16)  was  2nd  order  in  the  radial  direction).  These  boundary  con¬ 
ditions  can  be  obtained  by  considering  the  continuity  of  traction  and  velocity  across  the  mandrel/ 
composite  interface 

6n(x  =  \,y)  =  E'mg(y)  ,  (21) 

with  E'm  as  the  effective  modulus  of  the  mandrel,  and  utilizing  the  equilibrium  of  the  last  layer 
keeping  the  hoop  stress  constant  and  equal  to  the  applied  tension 


gJ*  =x„,y)  =  - 


T0dx0 
x0  dy 


with  T0  as  the  non-dimensional  applied  tension. 


(22) 


By  utilizing  (12)  and  (18),  the  boundary  condition  at  x  =x„  becomes 
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T0  dx0 


_dxa 

w,(x«,y) 

i  giy) 

j.  r/A  ^ 

\dwl+nwA 

~  dy 

L  cm 

*oGm_ 

K  dx  x  J 

(23j 


while  the  boundary  condition  at  x  =  1  becomes 


T0dx0  r’ 

x0  dy  J i 


m)~Tat; — p — *■' 


di, 


(24) 


The  fact  that  the  variation  of  the  outer  radius  with  time  is  unknown  provides  an  additional 
difficulty  in  wet  winding  analysis  which  is  absent  in  the  analysis  of  either  linear-elastic  or  viscoelastic 
dry  winding.  This  difficulty  is  related  to  the  fact  that  in  wet  winding  satisfaction  of  the  boundary 
conditions  leads  to  an  eigenvalue  problem  involving  the  outer  radius.  This  dependence  implies 
that,  unlike  the  linear  dry  winding  problems,  superposition  of  solutions  in  wet  winding  corre¬ 
sponding  to  different  outer  radii  leads  to  an  error  in  satisfying  force  equilibrium. 

3.  Numerical  Formulation 


One  major  source  of  biased  errors  during  integration  of  ( 1 6)  in  time  is  the  lack  of  equilibrium 
at  the  end  of  each  increment.  In  order  to  eliminate  this  type  of  error,  ( 1 5)  is  modified  in  the  following 
manner 


dorr(x ,  y  +  Ay )  (Jjx  ,y  + Ay)-  am(x  ,y  +  Ay) 
dx  +  x  + 


Ay 


dorr(x,y)  orr(x,y)-om(x,y) 
dx  +  x 


=  0 


(25) 


with  the  stress  rates  at  time  y  +  Ay  being  related  to  radial  velocity  at  that  time  through  the  constitutive 
relations  (12,13).  The  term  inside  the  bracket  represents  the  residual  body  force  at  time  y  +Ay  due 
to  the  lack  of  equilibrium  at  time  y.  The  above  equation  can  be  also  interpreted  as  the  force 
equilibrium  at  time  y  +  Ay  in  terms  of  the  stress  rates  at  that  time  and  the  total  stresses  at  the  previous 
time. 


The  form  of  the  modified  rate  equilibrium  equation  (25)  readily  renders  itself  to  Euler 
backward  method  (an  implicit  algorithm)  for  integrating  in  time.  A  finite  difference  scheme  was 
then  used  to  replace  the  spatial  derivatives.  Since  the  solution  of  the  problem  requires  the  ability 
to  follow  the  history  of  winding  from  beginning^  =0)  to  end  y  =  yf,  the  size  of  the  finite  difference 
mesh  was  sequentially  increased  as  time  increased,  see  Fig. 2.  Due  to  the  change  in  the  geometry 
which  is  induced  by  the  resin  filtration,  the  sequential  increase  in  mesh  size  resulted  in  non-uniform 
mesh  spacing.  By  applying  Richardson  extrapolation  (Dahlquist  and  Bjorck;  1974)  to  the  midpoint 
rule,  the  first  and  second  order  spatial  derivatives  of  the  interior  nodes  can  be  obtained  as 
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a_ 

dx 


~  dm(  )m  +  i  !  fn(  )m  -  1  +  Sm(  )m  , 


(26) 


I  l^m-l(  )m  +  l  +  ^*(  )m-l  (^m  fyn  -  1 )  (  )»J 


where 


4  = 


K-x 


»  /«  =  — 


«*=■ 


I  *  m  j  *  ~m  I 

‘m  ‘m 

hm=*m  +  l-Xm  , 


(27) 


(28) 


Similar  expressions  can  also  be  derived  for  the  boundary  nodes.  Due  to  the  large  variations 
that  the  function  F(<J>,)  can  experience,  me  governing  equation  for  w,  will  be  stiff.  In  order  to  reduce 
the  limitation  that  this  stiffness  imposes  on  the  time  increment  for  stable  integration,  we  proceed 
in  the  following  manner.  First,  we  integrate  the  evolution  equation  for  the  fiber  volume  fraction 
(11)  between  the  time  y  >  yn  and  y  =  y„.  The  result  will  be 


In 


rh(.x,y) 


-r 

>  >  *>« 


dws  w/| 

Lar+rr 


(29) 


By  utilizing  the  above  relation,  one  can  then  determine  the  following  relation  for  the  spatial 
derivative  of  fiber  volume  fraction 


0<t>,  (.*,}’) 
dx 


h(x,y)d$,(.x,yK) 
4>,Cx,yJ  dx 


4>. 


(  D2 


A 


d  Wj  1  dws  Wj 

^TT  +  -^ - 2  dy 

dx2  X  dx  X 2 ) 


(30) 


The  final  step  in  obtaining  the  finite  difference  eq  uation  of  the  original  PDE  will  be  to  substitute 
(30)  into  ( 1 6),  and  then  into  (25).  By  utilizing  the  midpoint  rule,  integration  of  the  resultant  equation 
between  the  time  y  =  yn  and  y  =  yn  + ,  leads  to 
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~ 

( 
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\ 
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V 
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1  rSn+1  gH.  1 
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,^(1  +  1  a.J  +  A  v  I^mCC(m  +  l)  +  em+  <Kr(m)  +  fnSCr(m-\) 


**Ay  G”"  Ay 
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,  for  2<m  <n  ,3<n 


(31) 


where 


2^-1  2/im  -2(/im  +  /im_,) 

=-7 —  +  <4  .  cm=  —  +fm  ,  6m  = - - - +  en:  . 

and  the  superscripts  specify  the  time  and  the  subscripts  specify  the  node  numbers. 
The  boundary  condition  (23)  is  replaced  by  the  following  difference  equation 


=!t7^|IV-T  -"*"*'1  +f.\*,W k::+<&  +-^-w;"+ w,"']  ,od 

TttI  ay  Cl)  U„+l  -*(l  +  l  *(1  +  1 


with 


d  = 


2hn  +  hn _ , 


+^»-l  7 

e  = — ~ —  ,  /= 


hn(hn  +  hn- 1)  ’  M.-1  ’  K  - \{K  +  K  - 1) 


The  finite  difference  equation  of  the  original  PDE  (31),  and  that  of  the  boundary  condition 
(32)  can  be  used  to  obtain  the  following  systems  of  equations  for  the  nodal  velocities  at  time 
>’„  +  i  =  v„  +  Ay 

AW”"  =B+Rg”"  ■  (33) 

where  W”"T  =  \W2 . VVJ"  +  1 

The  coefficients  of  the  tridiagonal  matrix  A ,  and  the  column  matrices  B  and  R  are  giver,  in 
the  Appendix.  Based  on  the  above  finite  difference  formulation  of  the  problem,  the  strategy  for 
the  computation  is  as  follows: 
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a)  For  the  first  time  increment,  by  assuming  linear  spatial  variation  of  the  velocity  of 
solid  the  boundary  conditions  (23)  and  (24)  are  used  to  determine  the  nodal  values 
of  ws. 


b)  During  the  subsequent  increments,  the  nodal  velocities  are  obtained  through  a 
predictor-corrector  method  with  each  pass  consisting  of: 
bl)  Using  (33)  to  obtain 

W**'  =A~'  [B+RgH  +  x]  ,  (34) 


b2)  Utilizing  the  boundary  condition  (24)  to  obtain  the  interface  velocity  gn  + 1. 
b3)  Determining  the  nodal  velocities  from  (34),  and  hence  the  increments  of 
stress,  pressure  and  fiber  volume  fraction. 

We  used  the  total  equilibrium  equation  in  the  radial  direction  to  obtain  the  following  estimate 
of  the  error  associated  with  the  numerical  calculation 


error  = 


0eeL 


(r°rr)  °09 


(35) 


4.  Discussion  and  Results 


In  this  section  by  assuming  specific  forms  for  the  functions  F(§s)  and  G  (<}),),  we  will  analyze 
response  of  the  purposed.  These  forms  are  given  as  follows 


(35) 


G(<t>J  =  G 


(36) 


with  6"  as  the  initial  fiber  volume  fraction  of  the  layer  which  is  being  wound  onto  the  mandrel, 
E°C°„  as  the  initial  stiffness  of  that  layer,  and  G°k°  as  its  permeability.  The  constant  F°  is  such 
that  F($a)C°rr  will  be  the  radial  stiffness  of  the  fiber  network  with  the  maximum  possible  volume 


fraction. 


The  form  for  G  (<j>J,  see  Fig.3,  is  chosen  based  on  the  expression  derived  by  Gebart  (1990) 

who  showed  that  the  variation  of  the  transverse  premeability  of  the  mixture  with  fiber  volume 
fraction  is  fundamentally  different  from  the  variation  of  the  axial  permeability.  The  form  of  the 
function  F{§s)  is  chosen  so  that  the  radial  resistance  will  initially  remain  almost  constant,  see  Fig.3. 
The  above  forms  have  the  interesting  feature  that  the  drop  in  the  permeability  is  faster  than  the  rise 
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in  the  radial  stiffness.  This  implies  that  based  on  the  assumed  forms  for  the  functions  F (<[>,)  and 
G  (<}),)  even  though  compaction  can  be  significantly  reduced  due  to  the  decrease  in  the  permeability, 
the  radial  stiffness  can  still  be  small. 

The  material  constants  and  the  value  of  the  parameters  associated  with  the  winding  process 
are  given  in  the  table  below.  In  order  to  reduce  the  amount  the  computation  involved,  we  have  only 
considered  winding  onto  a  rigid  mandrel. 

Initially,  the  variation  of  the  radial  stiffness  with  fiber  volume  fraction  was  neglected,  see 
figure  4-8,  however  the  permeability  of  the  mixture  was  allowed  to  vary  with  <(>f.  This  was  done 
so  as  to  investigate  the  influence  of  the  radial  stiffness  on  fiber  compaction.  The  time  increment 
for  each  calculation  was  At  =  0.2  sec.  Time  increments  greater  than  this  lead  to  significant  error 
accumulation  and  eventual  instability  of  the  result.  The  results  in  table  2  present  the  maximum 
error  in  equilibrium  which  happen  to  occui  in  the  outer  node. 

The  results  in  figure  4  shows  the  variation  of  the  fiber  volume  fraction  through  the  thickness. 
As  expected,  this  figure  indicates  that  increasing  the  radial  stiffness  of  the  fiber  network  suppresses 
fiber  compaction.  The  interesting  feature  of  this  figure  is  that  for  cases  where  the  initial  resistance 
of  the  fiber  network  is  sufficiently  small  so  as  to  allow  significant  compaction,  there  are  two 
boundary-layers:  one  close  to  the  mandrel-composite  interface  and  the  other  also  close  to  the  outer 
layer.  In  between  these  two  layers  fiber  volume  fraction  is  nearly  constant4.  Figure  5  presents  the 
details  inside  the  boundary-layer  close  tox  =  1.  As  the  results  indicate,  the  boundary-layer  becomes 
more  diffused  as  the  radial  stiffness  increases.  In  interpreting  this  figure,  it  is  important  to  note  that 
the  length  of  the  boundary-layer  is  on  the  order  of  the  thickness  of  an  individual  composite  layer. 
Thus  keeping  in  mind  that  fiber  diameter  is  on  the  order  of  10pm,  adopting  the  continuum  approach 
in  the  manner  presented  in  this  study  should  be  viewed  with  caution.  Apart  from  the  inner 
boundary-layer  which  is  dependent  on  the  stiffness  of  the  mandrel,  the  results  presented  in  figures 
4  and  5  imply  that  the  profile  of  the  fiber  volume  fraction  can  be  separated  into  to  a  passive  and  an 
active  zone  of  consolidation.  Inside  the  passive  zone  the  compaction  due  to  the  winding  tension 
has  reached  its  saturation  limit,  while  inside  the  active  zone  the  compaction  is  still  taking  place. 

Figure  6  show  the  variation  of  the  hoop  stress  in  the  composite  cylinder.  The  results  indicate 
that  the  drop  in  the  hoop  stress  (loss  in  tension)  increases  by  increasing  the  thickness  of  the  cylinder. 
The  loss  in  tension  is  more  dramatic  when  the  radial  stiffness  of  the  fiber  is  small.  The  drop  in 
hoop  stress  can  be  explained  by  observing  that  each  time  a  layer  is  added  to  the  cylinder  there  will 


4  Numerical  results  indicate  a  slight  linear  variation  of  the  fiber  volume  fraction  with  radial  posi¬ 
tion.  The  slope  of  this  inner  layer  increases  as  the  radial  resistance  of  the  fiber  increases. 
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be  a  compressive  hoop  strain  increment  and  a  compressive  radial  strain  increment.  These  com¬ 
pressive  strain  increments  result  in  compressive  hoop  and  radial  stress  increments.  The  larger  the 
magnitude  of  anisotropy  ( the  larger  X),  the  more  compressive  the  hoop  stress  increment.  Thus  as 
the  degree  of  anisotropy  increases,  the  loss  of  the  initial  winding  tension  increases.  The  results  in 
Fig.  6  also  indicate  that  in  the  case  where  the  radial  stiffness  is  very  small  (case  1),  the  hoop  stress 
w  ill  be  almost  zero  in  the  interior  of  the  cylinder. 

Figures  7  and  8,  respectively,  present  the  variation  of  the  radial  stress  and  the  pressure  through 
the  thickness  of  the  cylinder.  The  results  in  general  indicate  that  increasing  the  thickness  of  the 
composite  leads  to  increases  in  the  level  of  the  compressive  radial  stress  and  the  resin  pressure. 
Although  the  pressure  distribution  only  shows  the  existence  of  outer  boundary  layer,  the  actual 
numerical  results  do  indicate  both  boundary  layers.  The  results  also  indicate  that  for  a  given  number 
of  layers,  increasing  the  radial  stiffness  leads  to  a  drop  in  the  resin  pressure  and  an  increase  in  the 
level  of  the  compressive  radial  stress.  When  the  radial  resistance  is  sufficiently  small  that  significant 
fiber  compaction  can  take  place,  increasing  the  thickness  of  the  cylinder  leads  to  the  expansion  of 
the  passive  zone  of  consolidation.  This  can  be  seen  from  figure  8  where  increasing  the  number  of 
layers  leads  to  the  enlargement  of  the  region  in  which  pressure  is  constant5. 

Figures  9  through  14  present  the  results  of  the  numerical  simulation  in  which  the  radial  stiffness 
was  allowed  to  vary  with  the  fiber  volume  fraction.  The  details  pertaining  to  the  limits  on  the  radial 
stiffness,  the  time  increment,  and  the  maximum  error  of  each  calculation  is  shown  in  table  3. 

Figure  9  indicates  that  inclusion  in  the  analysis  of  the  variation  of  the  radial  stiffness  with  the 
fiber  volume  fraction  leads  to  an  increase  in  the  length  of  the  inner  boundary  layer.  Furthermore, 
as  figure  15  indicates,  considering  such  an  effect  causes  a  non-monotonic  dependence  of  the  resin 
pressure  on  the  attainable  radial  stiffness.  This  is  such  that  increasing  the  maximum  radial  stiffness 
from  1.476  MPa  to  147.6  MPa  will  actually  increase  the  pressure  inside  the  passive  zone  of  the 
consolidation.  Aside  from  these  differences,  inclusion  of  the  variation  of  the  instantaneous  elasticity 
of  the  fiber  network  with  fiber  compaction  does  not  alter  any  of  the  other  characteristics  of  the 
solution  with  constant  radial  stiffness. 


5  Note  that  based  on  the  Darcy’s  law  (14)  when  the  pressure  is  constant,  the  filtration  flux  drops 
to  zero  and  hence  the  consolidation  stops 
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Conclusion 

By  utilizing  the  rate  form  of  the  continuum  mixture  theory,  the  wet  winding  of  circumfer¬ 
entially  reinforced  composite  cylinders  is  studied.  The  purposed  formulation  contains  two  important 
characteristics  of  this  manufacturing  process,  namely  the  variation  of  the  instantaneous  permeability 
and  elasticity  of  the  mixture  with  fiber  compaction.  For  the  case  of  winding  onto  a  rigid  mandrel, 
the  numerical  results,  which  are  obtained  based  on  a  finite  difference  approximation,  indicate  that 
when  variation  of  the  permeability  with  the  fiber  volume  fraction  is  taken  into  account  there  will 
be  two  distinct  zones  of  consolidation.  One  in  which  fiber  compaction  has  ceased  to  take  place  and 
the  other  inside  which  the  fiber  bundles  are  continuing  to  compact  under  the  influence  of  the  applied 
tension.  Numerical  results  further  indicate  that  the  existence  of  these  two  regions  does  not  depend 
on  whether  or  not  the  increase  in  the  radial  stiffness  of  fiber  network  is  taken  into  account. 

Based  on  the  present  results,  an  undesirable  feature  of  this  manufacturing  process  is  the  loss 
of  tension  in  the  passive  zone  of  consolidation  away  from  the  mandrel.  However,  the  results  indicate 
that  this  drop  in  tension  (hoop  stress)  decreases  as  the  radial  stiffness  increases.  Thus,  if  one  is  able 
to  increase  the  radial  stiffness  through  the  frontal  curing  of  the  layer  close  to  the  mandrel,  one  will 
be  able  to  reduce  the  drop  in  tension. 

The  results  for  the  cases  when  the  radial  stiffness  was  allowed  to  increase  from  a  relatively 
small  value  to  a  large  value  show  that  one  can  achieve  substantial  resin  pressure.  This  is  a  desirable 
outcome  since  it  prevents  void  formation  during  the  curing  process. 

Further  improvements  on  the  present  analysis  can  be  achieved  by  developing  more  effective 
algorithms  for  taking  into  account  singularity  of  the  velocity  equilibrium  equation.  This  enables 
one  to  analyze  the  winding  problem  for  arbitrary  evolution  of  the  elasticity  of  the  network. 
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Appendix 


The  elements  of  the  tridiagonal  matrix  A ,  and  the  column  vectors  B_  and  R  are  given  as  follows 
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The  above  expressions  are  valid  for  all  the  nodal  positions  except  the  boundary  nodes  where  one 
has  to  take  account  of  the  boundary  conditions.  Utilizing  (32),  the  equivalent  expressions  for  the 
nodes  m  =  2  and  m  =  n  are 
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Velocity  of  deposition  w, 


0.256  mis 


thickness  of  the  last  layer  H 

0.25  mm 

applied  tension  T0  (dimensional  value) 

10  [MPa] 

Outside  radius  of  the  mandrei  /?, 

0.5  m 

Hoop  Stiffness  of  the  fiber  C& 

222  [GPa] 

Ce/C 

0.204 

initial  fiber  volume  fraction  <S?° 

0.5 

Table  1-  The  constants  used  in  the  numerical  solution 


Radial  stiffness  C" 

Error  j 

147.5  Pa 

10.56  % 

1.475  [kPa] 

1.78  % 

14.75  [kPa] 

0.217  % 

1.475  [MPa] 

5.02  x  104 

0.1475  [GPa] 

-  4.52  x  10’5 

Table  2  -  The  maximum  error  for  the  cases  where  the  radial  stiffness  was  considered  constant 


-19- 


Table  3 


Initial  radial  stiffness  C°(4£) 

Final  radial  stiffness  C“(< |>a) 

A t  sec 

Error 

1.475  kPa 

1.475  MPa 

0.2 

1.789  % 

1.475  kPa 

14.75  MPa 

0.2 

1.80  % 

1.475  kPa 

147.5  MPa 

0.2 

2.46  % 

1.475  kPa 

1.475  GPa 

0.1 

2.23  % 
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Fig.2  Schematic  of  the  finite  difference  mesh 
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Fig.6  radial  variation  of  the  non-dimensional  hoop  stress  for  the  following  cases:  1) 
Cl  =  14.75  kPa  and  n  ~  100  layers,  2)  C°  =  1.475  MPa  and  n  -  50  layers,  3) 

Cl  =  1 .475  MPa  and  n  =  100  layers,  4)  C°  =  0.1475  GPa  and  n  =  50  layers,  5) 

Cl  =  0.1475  GPa  and  n  =  100  layers. 
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§  =  (R  -  R,)/^-  Rj) 


Fig.7  radial  variation  of  the  non-dimensional  for  the  following  cases:  1) 

Cl  =  14.75  kPa  and  n  =  100  layers,  2)  C"  =  1.475  MPa  and  n  =  50  layers,  3) 
Cl  =  1 .475  MPa  and  n  =  100  layers,  4)  Cl  =  0. 1475  GPa  and  n  =  50  layers,  5) 
Cl  =  0.1475  GPa  and  n  =  100  layers. 
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s=(R-Rj)/(Ft,-Ri) 

Fig.8  radial  variation  of  the  non-dimensional  pressure  for  the  following  cases:  1) 
C®  =  14.75  kPa  and  n  =  100  layers,  2)  C®  =  1.475  MPa  and  n  =  50  layers,  3) 

C°  =  1.475  MPa  and  n  =  100  layers,  4)  C®  =0.1475  GPa  and  n  -  50  layers,  5) 
C®  =  0.1475  GPa  and  n  =  100  layers. 
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£  =  ( R  -  R.)/(  R0-  R.) 

Fig.9  Profile  of  fiber  volume  fraction  within  the  inner  boundary-layer  for  the  follow¬ 
ing  cases:  1)  C*  =  1.478  =  const  kPa,  2)  1.478  kPa  <  C„  <  1.476MPa,  3) 

1.478  kPa  <C„<  14.76MPa,  4)  1.478  kPa  <  C„  <  147.6MPa,  5) 

1.478  kPa  <  C„  ^  1.476GPa.  The  number  of  layers  n  is  held  constant  at  100. 
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s  =  ( R  -  R.)/(  R0-  Rj) 

Fig.  10  Radial  distibution  of  fiber  volume  fraction  for  the  following  cases:  1) 
C“  =  1.478  =  const  kPa,  2)  1.478  kPa  <  C„  <  1.476MPa,  3) 

1.478  kPa  <C„<  14.76MPa,  4)  1.478  kPa  <  <  147.6MPa,  5) 

1 .478  kPa  ^C„<  1.476GPa.  The  number  of  layers  is  held  constant  at  100. 
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Fig.l  1  Profile  of  fiber  volume  fraction  within  the  outer  boundary-layer  for  the  fol 
lowing  cases:  1)  C"  =  1.478  =  const  kPa,  2)  1.478  kPa  £  C„.  <  1.476MPa,  3) 

1.478  kPa  <C„<  14.76MPa,  4)  1.478  kPa  <  C„  £  147.6MPa,  5) 

1.478  kPa  <C„<  1.476GPa.  The  number  of  layers  n  is  held  constant  at  100. 
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Fig.  12  Radial  variation  of  the  hoop  stress  for  the  following  cases:  1) 

C*  =  1-478  =  const  kPa,  2)  1.478  kPa  <C„<  1.476MPa,  3) 

1.478  kPa  <C„<  14.76MPa,  4)  1.478  kPa  <  C„  <  147.6MPa,  5) 

1.478  kPa  ^C„<  1.476GPa.  The  number  of  layers  is  held  constant  at  100. 
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Fig.  13  Radial  variation  of  for  the  following  cases:  1)  C°  =  1.478  =  const  kPa,  2) 
1.478  kPa  <C„<  1.476MPa,  3)  1.478  kPa  C„  <  14.76MPa,  4) 

1 .478  kPa  <  C„.  <  147.6MPa,  5)  1.478  kPa  <  C„  <  1.476GPa.  The  number  of  layers 
is  held  constant  at  100. 
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Fig.  14  Radial  variation  of  the  resin  pressure  for  the  following  cases:  1) 

C°„  =  1.478  =  const  kPa,  2)  1.478  kPa  $  C„  <  1.476MPa,  3) 

1.478  kPa  <C„<  14.76MPa,  4)  1.478  kPa  <  C„  <  147.6MPa,  5) 

1 .478  kPa  <  C,,.  <  1.476GPa.  The  number  of  layers  is  held  constant  at  100. 
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Fig.  14  Radial  variation  of  the  resin  pressure  for  the  following  cases:  1) 

C  =  1.478  =  const  kPa,  2)  1.478  kPa  <*C„<  1.476MPa,  3) 

1.478  kPa  <C„<  14.76MPa,  4)  1.478  kPa  <  Cn  <  147.6MPa,  5) 

1 .478  kPa  <C„<  1.476GPa.  The  number  of  layers  is  held  constant  at  100. 


